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Abstract 

The problem of shock induced ignition by a piston is addressed in the framework of Fickett's model for reactive 
compressible flows, i.e., the reactive form of Burgers' equation. An induction-reaction two-step chain-branching 
model is used to study the coupling between the energy release and the compressible hydrodynamics occurring during 
the shock ignition transient leading to a detonation. Owing to the model's simplicity, the ignition and acceleration 
mechanism is explained using the two families of characteristics admitted by the model. The energy release along the 
particle paths provides the amplification of forward-travelling pressure waves. These waves pre-compress the medium 
in the induction layer ahead of the reaction zone, therefore changing the induction delays of successive particles. The 
variation of the induction delay provides the modulation of the amplification of the forward travelling pressure waves 
by controlling the residence time of the pressure waves in the reaction zone. A closed form analytical solution is 
obtained by the method of characteristics and high activation energy asymptotics. The acceleration of the reaction 
zone was found to be proportional to the product of the activation energy, the ratio of the induction to reaction time and 
the heat release. This finding provides a theoretical justification for the previous use of this non-dimensional number 
to characterize the ignition regimes observed experimentally in detonations and shock induced ignition phenomena. 
Numerical simulations are presented and analysed. Both subsonic and supersonic internal flame propagation regimes 
are observed, consistent with experiment and previous reactive Euler models. 

Keywords: Shock induced ignition, Fickett model, method of characteristics, activation energy asymptotics, 
amplification mechanism, coherent pressure wave amplification, Clarke equations 



1. Introduction 

The present study focuses on the problem of shock induced ignition by a piston suddenly accelerated into a 
reactive medium. This problem has attracted significant interest in recent years due to its universality in capturing 
the canonical processes occurring in reactive media in the presence of gas-dynamic disturbances, i.e. pressure waves 
[Q] [2l[3l ffl [6l [71 [U[9l [TO). The problem is also equivalent to shock induced ignition by a reflected shock ATI . 

Of particular interest to us is to determine the parameter that controls the different one-dimensional ignition 
regimes that can be realized experimentally ifTTll . For example, Meyer and Oppenheim first suggested that the sen- 
sitivity of the induction time to temperature fluctuations, with respect to the duration of the exothermic power pulse 
controls the ignition process [12 ] • Mathematically, this can be written as the product of the reduced activation energy 
controlling the induction reactions and the ratio of induction to reaction times fT3lL namely, 

Und E a 

x= 't — W (1) 
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The same parameter was also identified in the stability analyses of self-propagating detonation waves (l4l[T5l[T6l and 
has also been used for correlating the ignition regimes and stability of cellular detonations [ 1 3 1 . Typically, a value 
of x ~ 10 corresponds to the onset of strong pressure waves driven during ignition behind the shocks and the loss of 
stability of detonation waves. However, the physical mechanism responsible for changes of the shock-ignition regime 
when this parameter is varied remain unclear. 
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To address the relative importance of the induction zone sensitivity and the duration of the exothermic stage, the 
present study assumes the medium's decomposition is via a two step chain branching model, with distinct induction 
and reaction stages. This is typical of most reactive gases and can be shown to rationally derive from a generic three 
step chain branching model ITTIIIH). The investigation thus builds upon the studies of Dold and Kapila (6), Sharpe 
ifTOl and Blythe et al. |4 ], who have addressed this problem in the past using the reactive Euler equations coupled with 
two and three- step chain-branching models. 

Although significant insight has been gained in these past studies, the present work focuses on the much simpler 
Fickett model for compressible reactive media (1911201 • This model is an extension of Burgers' equation, a prototypical 
model for inert compressible hydrodynamics, to a reactive medium. Due to its simplicity, this model can serve to 
study the basic mechanism without the complications of the highly non-linear reactive Euler equations. For example, 
this simple model has allowed us in the past to clarify the instability mechanism of pulsating detonations and has 
reproduced the universal period doubling route to chaos in travelling detonations l2T1l . Using this model, we attempt 
to determine in the simplest form possible how the^f parameter controls the dynamics of the ignition process following 
shock compression. 

2. Fickett 's reaction model 

Fickett's analogue model to one-dimensional inviscid reactive compressible flow is given by l20ll2Tll : 



The model is formulated in Lagrangian coordinates, x denotes a material coordinate and t represents time lfT9l[20l . 
The variable p has the meaning of density in the model. The term p in ^ has the meaning of pressure, see l20l . For 
simplicity, we choose the generic equation of state proposed by Fickett, that is 



where A r is a reaction progress variable ranging from (unreacted) to 1 (reacted). The model has the property that 
pressure increases with increasing density, or release of energy at constant density. Note that one recovers the inviscid 
Burgers' equation by setting Q = 0. A similar simplified reactive flow model was also proposed by Majda |22], with 
slight differences in formulation. 

In the present work, we adopt the two step reaction model proposed by Radulescu and Tang [ 21 ] in order to 
capture the main features of chain-branching reactions |23|. Following the leading shock, we assume the existence 
of a thermally neutral induction zone, whose duration depends on the medium's local compression. Following the 
induction stage, we assume an exothermic reaction that proceeds at a state-independent constant rate. The resulting 
generic induction-reaction model we propose is thus: 



The variables have been non-dimensionalized by the characteristic state behind the initial shock. Time is thus 
non-dimensionalized by the induction delay behind the initial shock. The parameter e is the inverse activation density 
of the induction reactions. The parameter K, through our non-dimensionalization, is the ratio of the induction and 
reaction times. The Heaviside function H(-) controls the timing of the onset of the second exothermic reaction, which 
starts when the induction variable A[ reaches 0. Ahead of the shock, Ai - 1 and A r = 0. We are also assuming that the 
induction reaction is activated by the passage of the inert leading shock. The system to be solved is thus ([2]), ^ and 
§5§ with the equation of state given by ([3}. 

For the piston initiation problem of interest, the piston path corresponds to x = (see Fig. [T]). In order to model 
a steady piston, we set p(x - 0, t) = 1. Without loss of generality, we assume the strong shock limit and take the 
quiescent gas ahead of the shock as p(x, t = 0) = 0, Ai(x, t - 0) = 1 and A r (x, t = 0) = 0. For a non-reactive 
medium, Burgers' equation with the prescribed initial and boundary conditions result in a constant strength shock, 



dtp + d x p = 



(2) 




(3) 



d t Ai = -H(Ai)e^) 
d t A r = r = K{l-H{A i )){\-A r Y 



(4) 
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Figure 1: Characteristic t - x plot of shock initiation for K=0.2, e=0.2; the red line represents the main shock front t s (x), the green line represents 
the fire trajectory t*(x); the solid lines are C + characteristics; the broken line is the analytical prediction for the fire trajectory given by {26) for 
t = 0(l). 



propagating with speed 1/2, followed by a constant state p = 1 [ 24 1. When the medium is exothermic, the energy 
which is deposited in the medium alters the medium's mechanical response; we wish to analyse this mechanical 
response below. 



3. Amplification mechanism 

We first solved the above system numerically by a finite volume technique described by Radulescu & Tang |2T1 . 
In the present study, a grid resolution of 1600 grid points per unit length was used. We also set Q = 1. 

Figure [T] shows an example of the solution obtained numerically for K = 0.2 and e = 0.2. The variation of p is 
shown in the accompanying Video 7, but can be deduced directly, as it is given by the slope of the pressure waves 
(see below). The trajectories of the leading shock (the red line in the figure), the onset of energy release, which we 
will henceforth call fire (green line), and the internal pressure waves (black lines) are shown in the t-x diagram. The 
pressure wave emitted at the onset of reactions at (x = 0, t = 1) marks the end of the inert solution. The subsequent 
solution consists of an accelerating shock and fire, which eventually form a detonation wave. 

The physical mechanism responsible for the amplification shown in Fig. [T]is best explained by first transforming 
the governing equations in characteristic form and integrating ^ from the time a particle crosses the shock at t = t s (x) 
to the time at which the exothermic reactions begin at the fire t = t*(x). We obtain: 



dp 1 _ 
— = -rO 

dt 2 ^ 
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By inspection of ([6]), we deduce that the pressure waves propagating forward with acoustic speed dx/dt = p 
amplify whenever they travel through a reactive field. This is the case of the pressure waves in Fig. [T] originating at 
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Figure 2: Space-time diagram illustrating the run-away process at early times. A C + pressure wave originating at the piston (x = 0) amplifies while 
travelling through the reaction zone (R). Its emergence into the induction zone modifies the post shock state from (S) to a (D). This shortens the 
ignition delay (t* - t s ). The shorter ignition delay in turn provides a more coherent amplification of C + pressure waves in the reaction zone. Zones 
O and P denote the initial and product states, respectively. 



the piston once the reactions have begun at t = 1. The rate of change of the pressure wave's amplitude varies with the 
heat addition Q at the rate r. Indeed, Figure [T] shows these characteristics accelerating forward. For a constant reaction 
rate, the pressure waves amplification is proportional to their residence time in the reaction zone. Once these waves 
exit the reaction zone, they propagate into the induction zone, compressing the unreacted material and strengthening 
the leading shock. 

The other key element in the amplification stage is the dependence of the induction delay on this pre-compression, 
via 0. For low e (i.e., high activation energy), the induction delay can be significantly reduced. The tendency towards 
lowering the ignition delay permits the pressure waves to remain in the reaction zone for longer times, as can be seen 
from Fig. ([T]). This provides a more coherent amplification of these pressure waves, which in turn enhances the 
compression in the induction zone and that of the shock. The feedback mechanism is thus between the pressure waves 
propagating forward and the modification of the onset of reaction, which modulates the amplification amount. This 
mechanism was also found central in controlling the detonation wave stability by Radulescu and Tang |2H . 

4. Asymptotic analysis 

The early transient dynamics illustrated in Fig. [T]and discussed above have been summarized in Figure [2] The 
arguments presented can now be posed mathematically in an asymptotic analysis, which assumes that the inverse 
activation energy is small (i.e., e « 1). The framework is similar to that used by Sharpe ifTOl for the reactive Euler 
equations. The analysis aims to capture the coupling between the amplification of forward travelling pressure waves 
and the resulting shortening of the ignition delay discussed above. We assume that the energy release rate is slow, i.e. 




(9) 
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This makes the solution accurate on time scales of order unity (i.e. induction time scales). The analysis is also valid 
for fast heat release, which can be obtained by rescaling the time variable for faster time scales ITOlL restricting the 
result to earlier times. In both cases, what we are interested in is the acceleration of the fire. We will also restrict our 
analysis to the early amplification occurring before the lead shock has been modified by the arrival of a pressure wave 
originating in the reaction zone. We expand the density and reaction progress variable in the form 

p(x, t) = pi(x, t) + ep 2 (x, t) + 0(e 2 ) (10) 
A r (x, t) = A rl (x, t) + eA r2 (x, t) + 0(e 2 ) (11) 
**(*) = t\(x) + 0(e) (12) 

Replacing fLO} , fTT] ) and fT2] ) in the governing equations and boundary conditions, we obtain to leading order: 

Pi (*,*) = 1 (13) 
A rl (x,t) = (14) 
t\(x) = 1 + 2x (15) 

This corresponds to the inert solution of a piston driven shock wave. 
With the 0(1) solution obtained, the 0(e) problem to be solved is: 

d l 0\ 1 dx 

jM + a 4)=^ q along Tt= l (16) 

dX r2 dx 

— along -=0 (17) 

p 2 (x = 0,0 = (18) 
A r2 (x,t<f) = (19) 

w*U) 

e p2 dt = 1 (20) 

Jt s (x) 

The resulting problem is one of linear acoustics, which is very similar in form to what is obtained for the reactive 
Euler case, i.e., the Clarke equations (3[T0|. According to ([16]), the forward facing pressure waves all propagate at the 
constant speed (of unity), but amplify at constant rate in the reaction zone due to energy release. The solution can be 
obtained by the method of characteristics. In the reaction zone, for t > t* (zone R in Figure [2]), the reaction progress 
variable can be immediately obtained. 

A r2 (x,t>f) = at-t*(x)) (21) 

Substituting ([21} into ( p~6] > and using the boundary condition ( [T8| ), we obtain the variation of density in the reaction 
zone: 

P2 = ^ (f\x) - 1) (22) 

Once the density is known in the reaction zone, the value of density in the disturbed region D (see Figure [2]) can be 
found by extending the C + characteristics into this region. In region D, A 2 = and the C + characteristic relation ( [T6| ) 
requires that density is constant along a characteristic. The constant is evaluated using ( [22] ) at the intersection of the 
C + characteristic with the path t = f(x). Denoting this reference position x re f(x, t) (see Figure[2]), we get: 



P2(x, i) = ^Y (t*(x ref (x, t) - l) 



(23) 



with x re f given implicitly from the trajectory of the C + characteristic, i.e. 

X X re j 



, r< . = 1 (24) 

t-t \X re f) 
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Region S is bounded by the shock t s (x) = 2x and the first C + disturbance originating at (x - 0, t = 1), given by 
t = l + x. In this region, using the method of characteristics, we get pi(x, t) = A^ix, t) = 0. The solution is now 
complete, and the integral relation ( [20] ) can now be written as 



-4(f(x ref (x,t)-l) dt _ x = (25) 



Jl+x 



This expression, along with ( [24] ), give the trajectory t*(x) implicitly. Adopting Sharpe's iterative strategy ifTOl to find 
t*(x), we substitute the first order approximation t\(x) in the integrand and solve for the correction appearing in the 
integral's upper bound. From the first iteration, we find the following approximation: 

**(*) * 1 + x + ^ log (1 + Q£x) (26) 

This expression is shown in Figure [T] it reproduces quite accurately the acceleration of the reaction path. Evaluating 
the acceleration at the origin, we obtain: 

This concludes our analysis, which shows quantitatively the feedback mechanism discussed in the previous section. 
More importantly, however, it shows explicitly that the acceleration of the reaction zone is given as the product of x 
and the heat release Q. 



5. Influence of \ an d ignition regimes 

To illustrate the effect of varying x on me process of ignition and detonation formation, we have conducted 
numerical simulations by varying both e and K. This permits us to determine the evolution of the flow field on longer 
time scales, when the analysis presented above breaks down (i.e. x » !)• Keeping e = 0.2 as above and increasing 
K to a value of 2 provides a more rapid amplification, as shown in Figure [3] The accompanying Video 2 shows the 
evolution of the p field. This corresponds tox = 10. As can be seen in the figure, a shock wave now forms inside the 
induction zone from the confluence of the forward travelling pressure waves. This is due to the stronger amplification 
of the pressure waves in the reaction zone. This subsequent motion of the internal shock and fire are now in phase, 
and propagate quasi- steadily. The model thus recovers the internal reaction wave dynamics observed experimentally 
rm and for the reactive Euler equations |10|, where both subsonic and supersonic waves can be established. The 
supersonic waves occur when the initial acceleration of the fire is sufficiently large. 

Further increase in the rate of energy release rate K leads to a more prompt amplification of the pressure waves. 
Figure|4]shows the evolution of the flow field when K = 5, i.e. x = 50. The accompanying Video 3 shows the evolution 
of the p field. The internal shock wave now forms at an earlier time. The stronger compression in the induction zone 
now gives rise to a significant reduction in the induction delay. The onset of energy release now becomes in phase 
with the internal shock motion, and the two slowly accelerate. The arrival of this internal shock to the lead shock 
transforms the lead shock into a detonation wave. 

This interior wave takes the form of a quasi-steady weak internal detonation. It is not self-sustaining, as the 
expiration of the induction delay essentially sets its trajectory. For reference, a Chapman- Jouguet internal wave 
propagating into the induction zone medium (p = 1) would have a speed of 2 |20li . Instead, we find that the speed 
of this internal detonation wave is approximately 1.5 when it encounters the leading shock. For even larger values of 
K, we see the same evolution, but the internal weak detonations take on velocities higher than the corresponding CJ 
values. 

We have further explored the effect of varying both e and K in order to determine if a unique value of x can be 
used to characterize the acceleration process. While this is so for^ = 0(1) through the analysis presented above, 
larger values of x yield internal shock waves, invalidating the analysis once the shock has formed. Figure [5] shows 
three different cases where both e and K are varied so as to keep a constant^ = 20. For all three cases, we see that the 
formation of the internal shock and subsequent amplification are very similar, in spite of large differences in reaction 
rates. This suggests that indeed the entire acceleration process can be characterized by the unique parameter^. The 
analysis of the acceleration of the fire in the presence of the internal shock is left for future study. 
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Figure 3: Characteristic t - x plot of shock initiation for K=2, e=0.2; same legend as with Fig. ^ additionally with blue line representing the end 
of the reaction layer. 
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Figure 4: Characteristic t - x plot of shock initiation for K=5, e=0.2; same legend as for Fig. [3] 
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Figure 5: Characteristic t - x plot of shock initiation for^ = 20, a) K = 2, 6=0.1. b) K = 4, e=0.2. c) ^ = 8, e=0.4. 



6. Concluding remarks 



Modelling the shock induced ignition using Fickett's model permitted us to clearly illustrate the mechanism of 
acceleration of the induced hot spot: the energy release along the particle paths provides an amplification of forward- 
travelling pressure waves. These waves pre-compress the induction layer ahead of the reaction zone, therefore chang- 
ing the induction delays. The variation of the induction delay then provides a modulation of the amplification of 
forward travelling pressure waves. The model simplicity also allowed us to obtain an analytical solution for the ac- 
celeration of the reaction zone, which was found proportional to the product of the activation energy, the ratio of the 
induction to reaction time and the heat release. This finding provides a theoretical justification for the previous use of 
this non-dimensional number to characterize the ignition regimes observed experimentally in detonations and shock 
induced ignition phenomena. 

Acknowledgements 

We wish to thank the financial support of NSERC through a Discovery Grant to MIR and the support of the 
H2CAN NSERC Strategic Network of Excellence. 

References 

[1] J. F. Clarke, R. S. Cant, Nonsteady gas dynamic effects in the induction domain behind a strong shock, Vol. 95 of Progress in Astronautics 

and Aeronautics, AIAA, New York, 1985, p. 142. 
[2] P. A. Blythe, D. G. Crighton, Shock-generated ignition - the induction zone, Proceedings of the Royal Society of London Series a- 

Mathematical Physical and Engineering Sciences 426 (1870) (1989) 189-209. 
[3] L. Bauwens, Ignition between a shock and a contact surface: Influence of the downstream temperature, Proceedings of the Combustion 

Institute 28 (2000) 653-661. 

[4] P. A. Blythe, A. K. Kapila, M. Short, Shock-induced chain-branched ignition, Proceedings of the Combustion Institute 32 (2009) 2371-2377. 

[5] J. F. Clarke, D. R. Kassoy, N. E. Meharzi, N. Riley, R. Vasantha, On the evolution of plane detonations, Proceedings of the Royal Society of 
London Series a-Mathematical Physical and Engineering Sciences 429 (1877) (1990) 259-283. 

[6] J. W. Dold, A. K. Kapila, Comparison between shock initiations of detonation using thermally sensitive and chain-branching chemical-models, 
Combustion and Flame 85 (1-2) (1991) 185-194. 

[7] T. L. Jackson, A. K. Kapila, Shock-induced thermal runaway, Siam Journal on Applied Mathematics 45 (1) (1985) 130-137. 

[8] G. J. Sharpe, M. Short, Ignition of thermally sensitive explosives between a contact surface and a shock, Physics of Fluids 19 (12). 

[9] G. J. Sharpe, M. Short, Shock-induced ignition of thermally sensitive explosives, Ima Journal of Applied Mathematics 69 (5) (2004) 493-520. 
[10] G. J. Sharpe, Shock-induced ignition for a two-step chain-branching kinetics model, Physics of Fluids 14 (12) (2002) 4372-4388. 
[11] R. A. Strehlow, Fundamentals of Combustion, International Textbook Company, Scranton, Pennsylvania, 1968. 

[12] J. Meyer, A. Oppenheim, On the shock-induced ignition of explosive gases, in: Thirteenth Symposium (International) on Combustion, The 

Combustion Institute, 1971, pp. 1153-1164. 
[13] M. I. Radulescu, The propagation and failure mechanism of gaseous detonations : experiments in porous-walled tubes, Ph.D. thesis, McGill 

University (2003). 

[14] M. Short, G. J. Sharpe, Pulsating instability of detonations with a two-step chain-branching reaction model: theory and numerics, Combustion 
Theory and Modelling 7 (2) (2003) 401-416. 

[15] H. D. Ng, M. I. Radulescu, A. J. Higgins, N. Nikiforakis, J. H. S. Lee, Numerical investigation of the instability for one-dimensional 
chapman-jouguet detonations with chain-branching kinetics, Combustion Theory and Modelling 9 (3) (2005) 385-401. 

[16] C. Leung, M. I. Radulescu, G. J. Sharpe, Characteristics analysis of the one-dimensional pulsating dynamics of chain-branching detonations, 
Physics of Fluids 22 (12) (2010) 126101. 

[17] P. A. Blythe, A. K. Kapila, M. Short, Homogeneous ignition for a three-step chain-branching reaction model, Journal of Engineering Mathe- 
matics 56 (2) (2006) 105-128. 

[18] G. J. Sharpe, N. Maflahi, Homogeneous explosion and shock initiation for a three-step chain-branching reaction model, Journal of Fluid 

Mechanics 566 (2006) 163-194. 
[19] W. Fickett, Detonation in miniature, American Journal of Physics 47 (12) (1979) 1050-1059. 
[20] W. Fickett, Introduction to Detonation Theory, University of California Press, Berkeley, 1985. 

[21] M. I. Radulescu, J. Tang, Nonlinear dynamics of self-sustained supersonic reaction waves: Fickett's detonation analogue, Physical Review 
Letters 107 (16). 

[22] A. Majda, A qualitative model for dynamic combustion, SIAM Journal on Applied Mathematics 41 (1) (1981) 70-93. 
[23] W. Fickett, W. C. Davis, Detonation: Theory and Experiment, Dover Publications, Mineola, N.Y., 2000. 
[24] G. B. Whitham, Linear and Nonlinear Waves, Wiley, New York, 1974. 



9 



